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Mechanisms of acoustic energy dissipation in heterogeneous solids attract much attention in view of 
their importance for material characterization, nondestructive testing, and geophysics. Due to the 
progress in measurement techniques in recent years it has been revealed that rocks can demonstrate 
extremely high strain sensitivity of seismo-acoustic loss. In particular, it has been found that strains 
of order 10 -8 produced by lunar and solar tides are capable to cause variations in the seismoacoustic 
decrement on the order of several percents. Some laboratory data (although obtained for higher 
frequencies) also indicate the presence of very high dissipative nonlinearity. Conventionally discussed 
dissipation mechanisms (thermoelastic loss in dry solids, Biot and squirt- type loss in fluid-saturated 
ones) do not suffice to interpret such data. Here, the dissipation at individual cracks is revised 
taking into account the influence of wavy asperities of their surfaces quite typical of real cracks, 
which can drastically change the values of the relaxation frequencies and can result in giant strain 
sensitivity of the dissipation without the necessity to assume the presence of unrealistically thin 
(and, therefore, unrealistically soft) cracks. In particular, these mechanisms suggest interpretation 
for observations of pronounced amplitude modulation of seismo-acoustic waves by tidal strains. 



PACS numbers: 43.25.Ba, 43.25.Dc, 43. 25. Ed 

I. INTRODUCTION 

In recent years, much attention is paid to the so-called 
mesoscopic nonlinear elasticity 1 - of solids containing such 
structural features as cracks, contacts, intergrain aggre- 
gates of dislocations etc. that are small in the scale of the 
elastic wave length. Quite often the quadratic nonlinear- 
elastic parameter /? for such materials can be 10 3 — 10 4 
in contrast to ft ~ 10° typical of ideal crystals of ho- 
mogeneous amorphous solids. The common feature of 
the above-mentioned structural features-defects is their 
very high relative softness compared with that of the sur- 
rounding homogeneous material. Thus the local strains 
at the defects are strongly (often by several orders of 
magnitude) enhanced. This fact results in strongly in- 
creased macroscopic elastic nonlinearity, which can be be 
instructively elucidated using rheological-level models^ 3 - 
that are in principle distributed and heterogeneous even 
if the properties of regions small in the wavelength scale 
are described. 

In many cases, adhesion or frictional effects are also 
localized at those soft defects, which makes the resultant 
nonlinearity hysteretio 1 -. In recent years increasing at- 
tention is also paid to nonlinear-dissipative properties of 
mesoscopic solids, which manifest themselves in rather 
pronounced variations in the dissipation of one elastic 
wave in the presence of another wave 4 ^ or an applied 
(quasi)static stress (although the very fact of pressure- 
dependence of dissipation or its dependence on the acous- 
tic wave amplitude in rocks has been known for years&£) . 
The most striking feature is that quite moderate strains 
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e ~ 10~ 6 — 1CP 5 arc able to cause variations in the decre- 
ment up to tens of percents or even several times^^ 1 ^, 
whereas the accompanying variations in the clastic mod- 
uli are on the order j3e and do not exceed 10 -2 — 10~ 3 . 

In addition to laboratory measurements, even more gi- 
ant strain sensitivity of the dissipation is indicated by 
some field data. For example, in experiments 1 ^- on long- 
range (357 km and 430 km) propagation of monochro- 
matic elastic wave produced by high-stability vibration 
sources operating at frequencies of 5 — 7 Hz, the accuracy 
of the measurements was sufficient to single out periodic 
variations in the received-signal parameters which were 
well-correlated with the periodicity of the lunar-solar 
tides. For tidal strains in the Earth-crust, the character- 
istic amplitude 12 is ~ 1CP 8 , whereas in observations^ 
the variations of the received-signal amplitude amounted 
to 2 — 4 percents and 1 — 2 degrees for the signal phase. 
These values correspond to the path-averaged relative 
variations in the elastic modulus AE/E ~ 10~ 5 and ab- 
solute variations in the decrement AO ~ 3 • 10~ 5 which 
are of the same order of magnitude. Assuming for the 
average value of the decrement ~ (0.3 — 1) • 10~ 2 (i.e., 
Q = tt/0 ~ 100 — 300), we estimate that the relative vari- 
ation in the decrement is AO/0 ~ 0.3 • 10 ~ 2 and is over 
two orders of magnitude greater than AE/E ~ 10~ 5 . 

Taking into account that the tidal strain £0 ~ 10~ 8 , 
we obtain the path-averaged estimate for the quadratic 
nonlinear parameter [3 = (AE/E)/eq ~ 500 — 700, 
which is not so high as /? ~ 10 5 reported, for exam- 
ple, in the pioneering observation s 13 ' 14 of the tidal varia- 
tions in the elastic-wave velocities. The above estimated 
smaller value of (3 is not surprising since in the long-range 
experiments^ 1 - the wave path reached depths of several 
tens of kilometers where the high pressure of the over- 
burden rock layers closed soft cracks and significantly 
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reduced the path-averaged nonlinearity of the rocks. 

Much higher (comparable with the dat a 13 ' 14 ) values 
of f3 were observed in cross-well experiments 1 ^ with a 
high-stability downhole seismo-acoustic source operating 
at a frequency of 167 Hz, the propagation distance 360 m 
with the estimated wave velocity along the path about 
3000 m/s. In worki£ the tide-induced variations for the 
wave phase were about 0.05 rad. and elastic-modulus 
variations AE/E ~ 10~ 3 that indicated the nonlinear- 
ity parameter /3 = (AE/E)/e ~ (1 - 2) ■ 10 5 like in 
work s 13 ' 14 . For the wave amplitude, its tide-induced vari- 
ations were about 10% and corresponded to the absolute 
variations in the decrement A9 ~ (2 — 5) • 10~ 3 , which 
means that the strain-sensitivity of the rocks was so gi- 
ant that the tidal strains 10 -8 were able to produce the 
relative variations in the decrement A9/9 ~ 10~ 2 — 10 _1 . 

Similar estimates for A8/9, although less directly, are 
confirmed by the field observations of the tidal modula- 
tion of the intensity of received endogenous seismic noises 
at several observation sites at the Kamchatka peninsula 
and in Japa n 16 ' 17 . Normally, the weak-amplitude noise 
with strains 10" 10 — 10~ 12 was recorded by a sensitive 
narrow-band receiver around a frequency of 30 Hz. Co- 
herent averaging (from several weeks to several months) 
of the noise envelope made it possible to reliably sin- 
gle out periods of individual solar and lunar-tide com- 
ponents in the noise-intensity modulation with a typical 
depth ranged from 2 — 3 to 6 — 8%. Taking into account 
that for rather weak tidal strains eo ~ 10~ 8 , their direct 
influence on the rock fracturing and the accompanying 
seismo-acoustic emission does not look very probable, the 
observed modulation can readily be explained^ by the 
tidal modulation of the effective size of the region from 
which the signal at the receiver is collected. This size is 
determined by the characteristic damping length for the 
noise. Thus, the relative variations in the intensity of 
the received noise should be proportional to A9/6. The 
above obtained estimate A9/9 ~ 10~ 2 — lO^ 1 based on 
the independent direct measurements^ well agrees with 
the depth of the tidal modulation of the nois o 16 ' 17 . 

Thus various experimental data require an explana- 
tion of the giant value of strain-sensitivity of dissipation 
in mesoscopic solids. Since those strain-induced varia- 
tions in the dissipation are observed for very small am- 
plitudes of probing acoustic waves, for which absolute 
displacements at the microstructural defects (cracks and 
contacts) fall into essentially sub-atomic range, the re- 
sponsible mechanism(s) should not involve an activation 
threshold (unlike hysteretic mechanisms of frictional 6 -^ 
or adhesion origin 2 ^). 

As discussed in works 4,21 , pronounced strain- 
dependent dissipation in mesoscopic solids should arise 
due to combined action of purely elastic nonlinearity of 
the soft defects and conventional linear (i.e., viscous-like) 
dissipation that is also localized at the same defects be- 
cause of the locally strongly enhanced strain rate. This 
mechanism does not require a finite threshold (unlike es- 
sentially super-atomic displacements required for activa- 
tion of adhesion/frictional phenomena), although it can 
act in parallel with hysteretic mechanisms. 

In what follows, we consider physical realizations of 




FIG. 1. Schematically shown crack with wavy roughness at 
the interface which results in the creation of elongated (strip- 
like) contacts or waists inside the crack that may act as a kind 
of valve for the fluid flow is the crack is liquid-saturated. 



this mechanism that are relevant to solids containing dry 
and fluid saturated cracks. In both cases, the dissipation 
is threshold-less in amplitude: of thermoelastic origin in 
dry and viscous in fluid-filled cracks. In both cases the 
key role is played by the same geometrical features of real 
cracks for which corrugated surfaces (having wavy asperi- 
ties) are typical rather than smooth nearly plane-parallel 
form often used in the dissipation models. The analy- 
sis will be performed in the style of physical argumenta- 
tion used by Landau and Lifshita^ 2 - in the discussion of 
thermoelastic loss in polycrystalline solids and in work 23 . 
Such an asymptotic approach gives clear representation 
of the physics of the discussed phenomena and ensures 
quantitative estimates with a reasonable accuracy com- 
parable with that for formally exact solutions obtained 
for idealized (and thus approximate) models like elliptical 
cracks, etc. 



II. GENERAL CONSEQUENCES OF WAVY 
ROUGHNESS OF SURFACES IN REAL CRACKS FOR 
STRAIN SENSITIVITY OF DISSIPATION 

We emphasize the fact that unlike often assumed near- 
parallel geometry, for interfaces of real cracks, wavy 
forms are quite typical, which is confirmed by crack im- 
ages obtained by various methods and is in agreement 
with known models of crack initiation. Such initially co- 
inciding surfaces often are not simply separated in the 
normal direction but also exhibit certain tangential dis- 
placement and create inside the crack elongated "waists" 
(either nearly contacting or already contacting) as shown 
in Fig.[T] Inside liquid-saturated cracks, the so-created 
narrow waist can act for fluid flows as a kind of valve 
near which pressure gradients, flow velocities and the 
corresponding viscous dissipation in the crack should be 
localized. For dry cracks, the thermo-elastic dissipation 
should also be strongly localized and enhanced at the in- 
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ner contacts. In the vicinity of such wavy asperities, the 
local separation (or interpenetration) h of the crack sur- 
faces is significantly smaller than the average opening h 
of the crack. In contrast, the absolute variation in the 
average opening h and in the local separation h of the 
crack surfaces are practically the same, Ah w Ah. Due 
to this fact, the variations of the contact pre-strain (or 
the narrow waist opening) can be h/h ^> 1 stronger per- 
turbed than the average opening which determines the 
loss at the whole crack. 

After finding the acoustic loss at one crack the macro- 
scopic logarithmic decrement 8 can be found as the ratio 



2W a 



(1) 



where the energy Wq is dissipated during one period in 
a unit volume and W ac = Ke 2 /2 is the acoustic energy 
density stored in the material with the modulus K. In the 
order-of-magnitude estimates, it is not critical to specify 
which particular modulus is taken. Besides, we also can 
approximately consider that the presence of cracks does 
not drastically change the material's elastic modulus K. 



III. THERMOELASTIC LOSS AT INNER CONTACTS IN 
DRY CRACKS 

Unlike homogeneous materials for which thermoelastic 
dissipation of elastic waves is often negligible, in solids 
with microstructure (e.g., polycrystalline^) thermoelas- 
tic dissipation significantly increases due to the presence 
of small (compared with the elastic- wave length) hetero- 
geneities that strongly increase temperature gradients. 
For crack-containing solids, there is another factor which 
additionally strongly enhances the thermo-elastic cou- 
pling: the stress- and strain concentration at the crack 
perimeter as considered in work 2 ^. For the thermoelastic 
loss at the entire crack with characteristic diameter L, 
the thermoelastic dissipation exhibits maximum in the 
vicinity of the relaxation frequency /j, « k/ (2irpCL 2 ), 
where p is the material density, C the specific heat per 
unit mass, and k the thermal difiusivity. For millimeter- 
size cracks in rocks, this maximum corresponds to fre- 
quencies in the range 10^ 3 ..10~ 1 Hz. The analysis^ was 
based on exact solutions for the stress-field distribution 
near two-dimensional cracks represented as narrow ellip- 
tical cavities. However, even without specifying details of 
a particular crack model and estimating temperature gra- 
dients determined by the crack size in the low-frequency 
limit and by the temperature-wave length in the high- 
freqeuncy limit, it is possible to evaluate the elastic en- 
ergy loss using the approach similar to that used by Lan- 
dau and Lifshitz 22 for polycrystalline solids. 

We start from the thermal difiusivity equation for tem- 
perature variations T with respect to the mean value Tq 



dT k — d e 



(2) 



where 7 = ptK j (pC) is the Gruneisen parameter of the 
thermoelastic coupling, p? is the thermal expansion co- 



efficient, K is the bulk modulus, e is the material dilata- 
tion which for the present approximate consideration can 
be identified with strain in the field of a compressional 
wave. The acoustic energy loss due to irreversible heat 
flows can be found from the integral 22 



dW k , . 

-w = -%jw dv - 



(3) 



Estimating the gradients in the crack from Eq. © and 
evaluating integral ([3]) we obtain the following approxi- 
mate expressions for the elastic energy dissipated by the 
crack during one oscillation period, which well agree with 
the asymptotic forms of the results presented in2£: 



W^lck « 2irujTo(p' T K'/ K )L b e\ for u < uj l 



pCL 2 ' 
(4) 



WfZk « 27rT ^L[-^] 1 / 2 LV, for u » u L , (5) 
pL pCoj 



TF, 



crack 



2irT„(p^K 2 /pC)L A e\ for u w u L , (6) 



where oj is the circular frequency and is the character- 
istic circular frequency of thermal relaxation determined 
by the characteristic diameter L of the crack. 

Since we need Eqs. ©-© only for comparison with 
similar equations for thermoelastic loss at inner con- 
tacts in cracks, here, we only briefly outline their deriva- 
tion. First, note that for sufficiently low frequencies, the 
crack size L is much smaller than both the elastic wave 
length and the thermal wave length X t h — [/^/(Cpu))] 1 ' 2 . 
Therefore, it is the scale L which determines the gra- 
dients of the temperature variations in the elastic-wave 
field, \/T ~ T/L. Thus in Eq. ©, we can estimate 
that AT - f/L 2 and df/dt ~ uif. The condition 
L <C Xth is equivalent to w <C wl = ^/(pCL 2 ) and thus 
df/dt < AT, so that the first term in Eq. © can be 
omitted. Since de/dt ~ we, Eq.© yields the estimate 
for the temeprature variation T ~ lojTqC 'pL 2 / k. Then 
we estimate integral j2]) taking into account that the tem- 
perature gradient T/L is localized in the crack vicinity 
within the characteristic volume L 3 . Then one obtains 
dW/dt ~ w 2 To7 2 C 2 p 2 T 5 e 2 . Finally, for the energy dissi- 
pated over one oscillation period 2tt/u> we recover Eq.©. 

In the high-frequency limit when ui 3> ul — K / {pCL 2 ), 
in contrast, the first term dT/dt becomes dominant in 
Eq. I|2]). so that the the temperature variations are adi- 
abatic: T w 7To£ and Xth << L. In this case, the dis- 
sipation is mainly localized in the near vicinity of the 
crack edge, in which fairly universal asymptotic behavior 
of strain 2 ^ has the form e cr ~ s{r / V)~ x l 2 . Here, ra- 
dius r is counted from the crack edge (and the angular 
factor in this order-of-magnitude estimate is neglected). 
When estimating integral © we have to consider two 
regions r > X t h and r < Xth- It the latter region, the 
formally infinite strain s cr ~ e{r / L)^ 1 / 2 and strain gra- 
dient at r — > does not cause divergence of the inte- 
gral, since the amount of the generated heat is finite 
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and is smeared within the region r < Xth- As a result, 
both subregions r > Xth and r < Xth give functionally 
identical and approximately equal contributions, so that 
dW/dt ~ K,T j 2 L 2 /Xth and for the loss over one period 
we recover Eq. ([5]). Finally, Eq. © can be found as the 
crossover between Eq. ^ and Eq. (O at ui = ujl- 

To estimate the expected relaxation frequency, we 
can take parameters k ~ 0.015 W/cm/K and p ~ 2.6 
g/cm 3 typical of quartz-like rocks and consider cracks 
with L ~ 1 mm, which yields ft, = w/(2ir) ~ 10 _1 Hz. 
It has been estimated is work o 23 i 24 that thermoelastic 
losses described by Eqs. (j4])-([6]) can account for the ob- 
served values of Q-factor in dry rocks, especially for low 
seismic frequencies. On the other hand, it is seen from 
Eq. ^ that the height of the relaxation maximum is pro- 
portional to the cube of the characteristic crack size L. 
Therefore, to explain the typically observed in dry rocks 
values Q ~ 300 — 1000 for ultrasonic frequencies and even 
for frequencies of the order of 10 2 — 10 3 Hz, it is neces- 
sary to assume very high densities of very small cracks 
(of micro-meter scale), which does not look realistic. 

In this context, the above mentioned narrow inner con- 
tacts have much higher relaxation frequencies. However, 
the question arises whether such contacts (whose vol- 
ume is much smaller than that for the entire crack) can 
dissipate an appreciable amount of the the elastic-wave 
energy? The expressions for thermoelastic energy loss 
at the inner contacts in cracks can be obtained much 
like Eqs. dJ])-© taking additionally into account the lo- 
cal concentration of strain at the contacts (se o 26 ' 27 ): 



W^ t » 2TTUjT (fi 2 rK 2 /K)l 2 L c L 2 e 2 , for u < 



pCl 21 
(7) 



W c H J t « (2n/u J )KTo(pTK/Cp) 2 L c (L/l) 2 e 2 , for u » wj, 

(8) 

« 2nT (p 2 T K 2 /pC)L c L 2 e 2 , for w a u,, (9) 

Here, L c in the length of the strip-like contact, / is 
its width and the other notations are the same as for 
Eqs. d3|-(|H]), but the relaxation frequency uji is deter- 
mined by the width 1 « L of the strip-like contact. 
It is seen from Eqs. ©-© that for frequencies much 
lower that the relaxation frequency u>i, the loss is grow- 
ing as a linear function of u>, and for frequencies much 
greater than u>i, the loss is inversely proportional to oj. 
The asymptotic law uj^ 1 is similar to the high-frequency 
asymptotic for spheroidal voids^ 4 - rather than flat cracks. 

The most striking conclusion, which is seen from 
Eqs. ([T])-©, is that for a strip- like contact with a length 
L c ~ L, the magnitude of loss in the vicinity of the max- 
imum located at uji ;g> uj is determined by the size of the 
whole crack (since L C L ~ L 3 ). Consequently the maxi- 
mum loss at such contact is of the same order as for the 
conventionally considered maximum for the entire crack 
(compare Eqs. ((6]) and (|9])). In contrast, the positions of 
the maxima uji ~ n/(pCl 2 ) and u>l ~ K /(pCL 2 ) on the 
frequency axis can differ orders of magnitude. Thus one 
crack of size L with a strip-like contact L c ~ L of width 
I <C L can produce near the relaxation peak w; the ther- 
moelastic dissipation comparable with the contribution 
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FIG. 2. Schematically shown relative positions and heights 
of the thermoelastic relaxation peaks for the crack of size L 
as a whole (curve 1), for the inner strip-like contact of width 
I <S L and length L c ~ L (curve 2), and a peak similar to 
curve 1, but for a small crack of size I (curve 1 '). It is assumed 
that L = 1 mm and l/L ~ 10~ 2 . 



of (L/l) 3 3> 1 small cracks with the diameter /. Taking 
again an example of quartz-like rock, this difference is 
illustrated in Fig. [2] for L/l = 10 2 , which means that one 
such inner contact produces the same dissipation as 10 6 
small cracks of size /. 

Noting that asymptotic Eqs. (JT])-© agree with the fre- 
quency dependence of loss for a classical relaxator and 
assuming that the density of such cracks with inner con- 
tacts equals no we obtain from Eq. §T§ the expression for 
the decrement in the entire frequency range: 



27rTop 2 nKL c L 2 uj/loi 



P C 



1 + 



r"n 



(10) 



To find the overall decrement determined by the con- 
tributions of contacts with different parameters L, L Cl 
and / we have to integrate Eq. dTU|) over the distribu- 
tion n(L, L c , I), for which it is reasonable to assume that 
the distribution over the width I should be essentially 
independent of the distribution over L and L c , so that 
n(L,L c ,l) = n(L, L c )n(l) and the averaging takes the 
form 

acont _ ZvTofj^K f L c L 2 uj/uji 



P C 



1 



^ 2 n(L, L c )n(l)dLdL c dl 
(11) 

We note that a wide distribution n(l) can strongly 
smoothen the frequency dependence of the dissipation re- 
sulting in nearly-constant Q-factor even without assum- 
ing very high densities of tiny cracks. 

It is seen from the structure of integral (fTTj) that the 
height of the relaxation peak and its location on the 
frequency axis are quite independent since they are de- 
termined by essentially different parameters (L and L c 
for the peak height and I for the position). Evidently, 
for sufficiently small variation in the average strains and 
stresses, the width of the contacts can already be sig- 
nificantly perturbed, whereas the length of the contact 
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FIG. 3. Variation in the shape 6(uj) of the normalized re- 
laxation peak for an inner contact (curves 1 and 2) due to 
15% variation in the contact width. The variation A8 in the 
decrement has opposite signs at the opposite sides from the 
maximum and peak-to-peak excursion about 15% (curve 3). 



and the average opening of the crack can be only slightly 
affected. This can result in the displacement of the re- 
laxation maximum on the frequency axis almost without 
affecting its height as illustrated in Fig. [3] for the case 
of identical contacts. To estimate how strongly the fre- 
quency u!i sy K/(pCl 2 ) of the maximum of the thermoe- 
lastic loss at such a contact can be perturbed by vari- 
ation in the mean strain in the material, let us use the 
presented in Landau and Lifshitz 22 solution for the width 
of the contact between two compressed aligned cylinders: 



(WDF RR' 
\ 3tt R + R> 



1/2 



(12) 



where D = (3/4)[(l - cr 2 )/E + (1 - cj' 2 )/E], a and a' are 
Poisson's coefficients of the materials of the cylinders, E 
and E' are the Young moduli of the cylinders, and F is 
the specific force per unit length of the contacts. In the 
case of identical materials with the same accuracy which 
corresponds to Ejqs.tfT])— (fTTJj) . we can consider that D s» 
3/ (2E) and the force F per unit length can be expressed 
as F = F c /L c , where F c is the total force applied to the 
contact. Then from Eq. (|12l) we obtain 



r 



(16F c R)/(irEL c ). 



(13) 



For the force F c applied to the inner contact in a crack, 
an approximate expression F c /L c rs eEL 2 j(L + L c ) can 
be obtained^!, where e is the mean strain in the material. 
Then taking into account that 1/2 < L/(L + L c ) < 1 one 
obtains from Eq. (|12[) the following approximate expres- 
sions: 



I 2 



it L - 



— — sRL ~ — sRL, 

L r 7T 



16 L _ 8 _ Auji _ A(l 2 ) RL 

This expression clearly elucidates the reason of the ex- 
tremely high sensitivity of I 2 (and, consequently, varia- 
tions in u>i) with respect to the variation Ae in the mean 



strain. Indeed, the wavy asperities at the crack interfaces 
often have radius comparable with the characteristic di- 
ameter of the entire crack (R ~ L), whereas the width 
of the inner strip-like contact often does not exceed the 
average opening h of the crack. Therefore, taking into 
account the above discussed characteristic values of the 
crack aspect ratios, we conclude that the factor RL/l 2 in 
Eq. (fl"4|) can easily be as great as 10 6 -10 8 . This means 
that even variations in the mean deformation ~ 10~ 8 
(which are conventionally considered unable to apprecia- 
bly influence acoustical parameters of solids) actually can 
be able to significantly change the width of inner contacts 
and affect the near-contact acoustic loss (see Fig. [3]). 

Estimating the overall sensitivity of the thermoelas- 
tic dissipation to the mean strain it is necessary to take 
into account that conventionally considered thermoelas- 
tic losses at the crack as a whole 2 ^ also contribute to 
the total decrement. We have, however, to recollect that 
the characteristic relaxation frequencies cjj « n/pCl 2 
and ujl ~ n/pCL 2 are related as (L/l) 2 and the high- 
frequency thermoelastic loss at the whole crack decreases 
as {uju} L y x l 2 (see Eq.©) as illustared in Fig.[U Then 
near the relaxation maximum uji we found that the ratio 
of the contributions of the loss 9i oc {uj — uj{) at the inner 
contacts to the conventionally considered "global" loss 
9 glob (u> = uji) at the cracks as a whole has the form 



Ogiob (w = u) L ) n cr 



1/2 



ric 
n, 



L 

'T 



(15) 



where n cr and n cr are the densities of cracks with and 
without inner contacts, respectively. Since L/l 3> 1, it 
is clear that even a small portion (e.g., a few percents) 
of cracks with narrow strip-like inner contacts can en- 
sure near u> = u>i the same contribution to the decrement 
as the conventional "global" mechanism of thermoelastic 
dissipation at whole cracks. Therefore, even quite a small 
portion of cracks with contacts can provide the above 
discussed extremely high sensitivity of the dissipation to 
small variations in the mean strain in the material. 

Even putting aside somewhat exotic manifestations of 
the giant stress-sensitivity of the dissipation and consid- 
ering only its mean value, it can be emphasized that in 
the audible and ultrasonic ranges typical of laboratory 
studies, inner contacts can ensure a fairly strong con- 
tribution to dissipation comparable with the estimates 
obtained in work s 23 ' 24 for the "global" thermoelastic loss 
at a crack as a whole for lower frequencies, which showed 
a reasonable agreement with the experimental data. 



IV. VISCOUS LOSS AT FLUID-SATURATED CRACKS 
WITH INNER STRIP-LIKE "WAISTS" 

Now let us consider viscous dissipation in cracks con- 
taining liquids. There is general agreement that along 
with the dissipation due to "global" fluid flows in pore 
channels (Biot's mechanism), an important role in the 
elastic-wave energy dissipation in rocks belongs to the 
local (or "squirt") flows inside crack o 28 ' 29 . Since the ge- 
ometry of pore channels for the Biot flows rather weakly 
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depends on the average strain and stress, it is the squirt- 
type dissipation in relatively soft narrow cracks which 
demonstrates much higher sensitivity to the variation in 
the mean strain in the material. In this section we focus 
on very important new features of the squirt mechanism 
in the case of cracks with wavy asperities of the interface. 

To estimate the viscous loss in narrow cracks 
with nearly parallel surfaces we will use the integral 
expression^ for the rate of kinetic energy dissipation in 
a flow of a viscous fluid between parallel solid planes: 



dW k . 



Of 



dv x 
dy 



dV 



(16) 



In integral (IT6"|) over the volume of the flow, Y axis is 
orthogonal to the interface, v x is the projection of the 
liquid-flow velocity on the X axis, and r\ is the fluid vis- 
cosity. The integration should be made over the volume 
of the flow. Using Eq. pl))) we first estimate the amount of 
energy Wq dissipated in a unit volume over one period. 
The form of Wq(uj) can again be understood by using 
asymptotic estimates very similar to those in the previ- 
ous section considering now the linearized Navier-Stokes 
equation for the flow velocity v x 



dv x 

at 



dp d 2 v x 
dx ^ dy 2 



(17) 



where p is pressure. In view of similar structures of 
Eq. (|17l) and Eq. @ , the roles of the diffusive term and 
the term with the time derivative can be estimated in a 
very similar way. 

In order to better delineate the place of our mecha- 
nism in the context of conventionally discussed ones we 
first recall how the well-known properties of the "global" 
Biot loss in fluid-filled channels can be derived from very 
simple arguments. For low acoustic frequencies, the role 
of the first inertial term in Eq. (|17l) is negligible, so that 
\dp/dx\ ?» \rjd 2 v x /dy 2 \, which corresponds to a fluid mo- 
tion as a parabolic Poiseuille flow in the channel. Thus 
the velocity gradient is determined by the entire channel 
thickness H, d 2 v x /dy 2 ~ v x /H 2 , and is frequency inde- 
pendent. The fluid velocity v x in the Poiseuille flow is in- 
duced by the gradient of the acoustic pressure p ac , which 
is determined by the acoustic wave length A w oc uj~ 1 , so 
that v x oc dp/dx ~ p ac /^uj oc w. Taking into account 
that integral (|16p is evaluated over a frequency- indepen- 
dent volume of the flow, the corresponding loss over one 
period 2ir/uj is Wq oc (2tt/uj) ■ uj 2 , so that 8(uj)lf oc uj. 

With increasing frequency to the first inertial term in 
Eq. (|17[) increases and becomes comparable with (and 
even greater than) the viscous term. This happens when 
Iwpu^l w \i]v x /H 2 \, i.e., for the characteristic frequency 



UJ,. 



r,/(pH 2 



(18) 



which is often expressed in a more indirect way via the 
ratio of porosity to permeability instead of the use of 
the characteristic width of the channels. The representa- 
tion of uj c in form (fT5)) is especially physically clear and 
means that the viscous wave length \ V i SC = (2r\j puj) 1 / 2 
becomes comparable with H. For u) 3> uj c the veloc- 
ity gradient becomes localized near the walls within the 



viscous layer X V i SC -C H, which determines the flow- 
velocity gradient, and we have (dv x /dy) 2 oc {uj 1 / 2 ) 2 oc uj. 
Next, since the last viscous term in Eq. (|17[) now can 
be neglected, we see that the velocity v x becomes pro- 
portional to the acoustical pressure amplitude p ac with- 
out any frequency-dependent factor, since dv x /dt oc ujv x 
and dp/dx ~ p ac /^u uj both are proportional to uj. 
Besides, in integral Eq. (fl6| the intergration should be 
made only within the viscous layer with the thickness 
(where the gradient is localized) rather 



A, 



oc UJ 



-1/2 



than over the entire volume of the fluid flow. The re- 
sultant frequency dependence of integral (| 10[) has the 
form cj _1 / 2 (cj 1 / 2 ) 2 = uj 1 / 2 . The latter estimate should 
be multiplied by the factor 2-k/uj to find the energy loss 
Wq 1F during one period, so that for ui 3> ui c , we obtain 
the asymptotic frequency dependence for the decrement: 
9{lo) hf oc (2tt/uj)uj 1 / 2 oc uj^ 1 / 2 . Thus we recovered 
the well-known asymptotic dependences of the decrement 
9{uj) due to the Biot flows in the pore channels. 

Let us now turn to narrow cracks which are much softer 
objects and can exhibit stronger stress-dependence. The 
above considered arguments based on the Poiseuille-flow 
approximation in the low-frequency limit remain very 
similar for flows in cracks and predict the same asymp- 
totic behavior of the decrement 9(uj) lf oc uj. As was 
mentioned in work—, the localization of the shear-flow 
gradients within a narrow viscous layer ~ \ visc for suffi- 
ciently high frequencies may also be applied to cracks, 
but the characteristic frequency is much higher, since 
the crack opening h is much smaller than the diameter 
H of the pore channels. However, in contrast to rigid 
pore channels, cracks are much softer objects^, with 
a characteristic own modulus aK <C K. For cracks, 
the aspect ratio a ~ h/L plays the role of their soft- 
ness parameter: due to their planar geometry cracks are 
roughly a^ 1 times more compliant than the surround- 
ing matrix with respect to both compression and shear. 
Thus for crack-like compliant pores, another characteris- 
tic frequency of viscous relaxation appears, which physi- 
cally corresponds to the condition that under sufficiently 
high-frequency oscillation, the viscous resistance to crack 
shearing becomes comparable with the elastic response 
of the crack to shearing. Namely, under tangential dis- 
placement Ax of the center of the crack, its shear elastic 
reaction apAx/h should be compared with the tangen- 
tial viscous stress r]UjAx/h. This yields a characteristic 
frequency 



uj d ~ ap/rj, 



(19) 



which was pointed out, for example, in work^i. Now we 
note that for geophysical materials, typically the values 
of the shear modulus p and compression modulus K are 
of the same order and the corresponding ujd is relatively 
high. For example, for viscosity r\ = 10~ 3 Pa - s typical of 
water and modulus fi ~ 10 10 N/m 2 typical of rocks and 
quite narrow cracks with a ~ 10~ 3 — 10~ 4 we obtain u>d ~ 
10 9 — 10 10 rad./s, which evidently can be relevant only 
to utrasonic laboratory studies, but not to the seismo- 
acoustic data discussed in the introduction. 

One more characteristic frequency can be obtained in 
a similar way, but considering the viscous resistance of 
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a crack-like defect to the compression in the direction 
normal to its plane. This viscous resistance should be 
compared with the elastic reaction also in the normal 
direction. To estimate the normal viscous resistance we 
have to supplement Eq. (fTTj) with the continuity equation 



dvx 
dx 



dvy 
dy 



= 



(20) 



since in the first approximation the liquid can be consid- 
ered incompressible. For the discussed nearly Poiseuille 
flow, dv x /dx ~ v x /L and dv y /dy ~ v y /L, so that from 
Eq. (|20l) we obtain v x ~ (L/h)v y . Since the variation 
Ah of the crack opening in the acoustic field corresponds 
to v y ~ wA/i, we find that v x ~ (L/h)ujAh. Next, in 
Eq. (|17|) . we can similarly estimate that dp/dx ~ p/L 
and rjd 2 v x /dy 2 ~ r/v x /h 2 . Neglecting the first inertial 
term in Eq. (|17l) for the Poiseuille flow and equalling the 
above estimated terms in the right-hand side of Eq. (fl7|) 
we obtain the estimate for the pressure p msc induced by 
the viscous flow of the liquid in the narrow gap: 



P 



VV x (L/h 2 ) ~ r]Uj{L/hf{Ah/h), 



(21) 



For the application to real cracks it is important to under- 
stand how critical is the parallelism of the crack surfaces 
(which for real cracks is rather a rule than an exception) . 
Assuming that the crack opening has the values h\ and 
hi at the opposite sides we conclude that, in the extreme 
case when at one side the h\ = 0, the maximal non- 
papallelism angle is of the order h/L ~ a, where h is the 
average crack opening. It is clear that for estimate (l2~T1) . 
this small non-parallelism is not important. On the other 
hand, one may argue that if the planes are not exactly 
parallel then the tangential viscous force acting along the 
inclined wall has a component normal to the crack plane. 
The viscous stress r msc acting along the inclined plane is 
readily estimated as r visc ~ r\v x jh ~ i](L/h)uiAh/h, so 
that its normal component is r^ ISC ~ r msc a ~ TjuiAh/h. 
Comparing the latter contribution and Eq. ([2~Tj) for the 
pressure due to the viscous flow, we see that p msc is 
1/a 2 >> 1 times greater than t% isc , so that actually 
the contribution of the eventual non-parallelism is not 
significant. Thus we can use Eq. (|21[) for the normal vis- 
cous reaction even for cracks with not perfectly parallel 
surfaces. 

Now we have another question: what is the elastic reac- 
tion that should be taken for comparison with the above 
found value (f2~Tj) for the viscous normal resistance. In lit- 
erature different variants can be found. One variant is 
by analogy with obtaining Eq. (fl9|) to compare Eq. (|2Tj) 
with the normal elastic reaction of the crack related to its 
effective modulus, that is, the quantity ~ aK . The lat- 
ter value corresponds to the well known rule of thumb 32 
that the effective softness relatively to the matrix ma- 
terial for a crack approximately equals its aspect ratio. 
This comparison gives us some characteristic frequency 
often discussed in literature^ 



faces by Ah (see Eq. (|2Tj) ) becomes so high that the near- 
crack region of the size L (where the elastic stress was 
initially nearly released) experiences compression by the 
value Ah* ~ p msc /(KL). For oj ~ u K , the value of Ah* 
becomes comparable with the initial approaching Ah and 
strongly compensates the latter. If the crack is "instanta- 
neously" compressed, initially the liquid remains almost 
"frozen" and then is redistributed during the character- 
istic time 1/u>k when the elastic compression Ah* of the 
near-crack region gradually releases. For sufficiently high 
frequency u> S> lok , the actually displaced volume of the 
liquid is equal to a small fraction wjf/u of the volume 
that would be quasistatically displaced, so that if the 
flow structure remains of the Poiseuille type, the amount 
of the dissipated energy during one period decreases as 
{uj/ujk) 1 - This observation will be used below to esti- 
mate the absorption in the relaxation maximum. 

It can also be pointed out that besides the relaxation 
frequency (|22p another relaxation frequency for cracks 
filled by compressible fluids is also discussed 2 ^. Namely, 
in view of the fact that modulus Kf of the filling fluid 
typically is at least an order of magnitude smaller than 
modulus K , it can be argued that the surrounding rock 
in the first approximation compresses the fluid as "ab- 
solutely rigid" body. Thus comparing the pressure of 
the compressed fluid KfAh/h with the above found nor- 
mal viscous reaction (|21[) of the fluid flow, we obtain 
another characteristic frequency uj r ~ K fa 2 jr\ which is 
determined by the compressibility of the fluid and in con- 
trast to Eq. (f2"2"l) does not depend on the modulus of the 
rock. a This new relaxation frequency differs from uik 
given by Eq. (l2"2"l) by a factor of Kf/(Ka). For narrow 
cracks with a <C 1 and typical liquids like water or oil, 
for which Kf/K ~ 10 _1 , the factor Kf/(Ka) is quite 
large, so that u r 3> U)k- I n such a case, long before 
the frequency ui r would be attained, the fluid flow should 
already be strongly damped for u) > u>k because of the 
finite rigidity of the near-crack region as discussed above. 
Therefore, the initial low-frequency asymptotic oc w can- 
not be extrapolated till the frequency Lo r 3> u>k- Con- 
sequently, the new (and potentially stronger) relaxation 
maximum near uj r simply will not be formed for typical 
liquids (like water and oil), for which Kf/(Ka) ^> 1. 
Only for sufficiently high-compressible (gaseous) fluids 
with Kf/(Ka) < 1 this mechanism should form the re- 
laxation maximum near oj r . But in this case the relax- 
ation frequency uj r should necessarily lie lower than u>k ■ 

Less formally, the physical meaning of this type of re- 
laxation is that a finite time is required for a compressible 
fluid squeezed inside a crack to form the flow which re- 
distributes the fluid end equalizes the internal pressure. 

In any case, the character of such flows can notice- 
ably be perturbed only if the mean opening of the entire 
crack is perturbed. This means that the mean perturb- 
ing strain should already be comparable with a like for 
the above-considered "global" thermoelastic loss. 



luk ~ c?Kjr\. 



(22) 



Less formally, it can be said that the viscous flow pres- 
sure p msc arisen due to approaching of the crack sur- 



Finally, another type of the characteristic frequency related to 
the fluid compressibility was considered in papeii22. For this 
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In what follows, we consider a modified form of the vis- 
cous loss in cracks with wavy roughness of the surface. 
This mechanism relates to the conventionally considered 
viscous loss at cracks much like the thermoelastic loss at 
the inner contacts relate to the thermoelastic loss at a 
crack as a whole. Even for "normal" liquids, the new 
strongly shifted to lower frequencies and very intense re- 
laxation maximum is related to the fluid compressibility 
rather than to compressibility of rock in the crack vicin- 
ity. Besides, which is especially interesting, the stress 
sensitivity of such a mechanism can be additionally en- 
hanced orders of magnitude and can readily explain the 
experimental data discussed in the introduction. 

Consider a crack with a "waist" created by wavy as- 
perities as shown in Fig. 2] The notations used are clear 
from the figure. Here, quantities p characterize the 
pressure in the respective cross sections of the crack. 
Assuming again the Poiseuille character of the flow for 
sufficiently low frequencies, the nearly plane fluid flow 
q = L z J v x (y)dy in the main part of the crack and in 
the narrower region of the waist should be equal: 



h 3 L z P 2 - Pi = h 3 L z P 3 - P 2 
12-q L ' \2r\ I 



(23) 



where the meaning of the pressure gradients (P 2 — P\)/L 
and (P 3 - Pi) ft is clear from Fig.H Then Eq. (J23J) yields 



P2-P1 
P 3 -P 2 



h 3 V 



(24) 



so that for sufficiently narrow waist in the crack, 
(h 3 /h 3 )(L/l) < 1, we conclude that (P 2 - Pi) < (P 3 - 
P 2 ) « AP tot , that is, the total drop AP tot of pressure is 
mainly localized in the region of the narrow waist. Now 
imagine an instantaneous slight compression of the crack 
under an external action after which the fluid compressed 
inside the crack flows through the narrowing into the 
outer region to equalize the pressure. This flow will expo- 
nentially decrease with a characteristic relaxation time r. 
Consequently the total volume q T of the leaked fluid dur- 
ing this relaxation process is given by the integral 



q T = q I exp{—t / r)dt = qr, 
'0 



(25) 



where q is the flow value corresponding to the initial ex- 
cess of pressure inside the crack over the pressure in the 
outer space. Since the rigidity of the outer pore space is 
much greater than that of the thin crack, the change in 
the volume of this rigid pore fraction is negligibly small 
compared with the strain of the crack. Consequently, the 
outer pore pressure remains close to its equilibrium value 



characteristic frequency, the elastic wave length becomes compa- 
rable with the length of the liquid drop inside the crack, so that 
the liquid drop cannot be already considered as a lumped object 
in the scale of the elastic wave length. If the liquid occupies 
almost all crack of diameter L, this frequency is given by the 
condition ui ~ (K f / p) 1 / 2 / L, which corresponds to fairly high 
ultrasonic frequencies and will not be further discussed. 




FIG. 4. Schematically shown crack with a narrow strip-like 
waist created by wavy asperities. The dominant part of the 
total difference in the pressure \Pi — Pi\ in the fluid flow is 
localized near the narrow waist, so that \P\ — Pi\/\P?, — P2I ~ 1. 
Pi is close to the equilibrium pressure Po in the outer pore 
channel. 



Po an d we can consider that P3 — P 2 
Thus the initial value of the flow 



P-P^Pi-Po. 



h 3 L z P 3 - P 2 h 3 L z P 1 - P Q 



12rj 



I 



12rj 



I 



(26) 



Since the variation \ i n the volume of the liquid and 
the pressure P are related via the bulk modulus Kf of the 
liquid, P = xKf, then the total liquid volume q T = qr 
displaced in order to equalize the pressure can be written 
via the initial volume Co of the fluid-filled crack as 



C (xi-Xo) = C Kf(P 1 -P ), 



(27) 



Equations ([26]) and ([27]) yield 



h 3 L z K f (xi 



Xo) 



12r] 



I 



Co(xi-xo). (28) 



Assuming that the sizes of cracks along X and Z direc- 
tions are of the same order L z ~ L, the crack volume can 
be estimated as Co f» hi? , so that from Eq. ([28]) we read- 
ily obtain the relaxation time r and its inverse value, i.e., 
the relaxation frequency uj r (tilde denotes that we con- 
sider the crack with a narrowing) : 




12r? 



h\ 3 I 



7 T- (29) 



o?K f \hJ L 



Here, we intentionally singled out the squared aspect ra- 
tio a 2 for convenience of comparison with the crack with- 
out waists^. The absence of the narrow waist can be 
interpreted as I ~ L and h ~ h, i.e., (h/h) 3 l/L ~ 1 so 
that Eq. ([2"TJ]) reduces to 



a 2 K 



f 



127] 



and r 



12r7 
a 2 K 



f 



(30) 



Within the accuracy of its derivation Eq. ([30]) coincides 
with the result r w 8r)/a 2 Kf obtained in work2& for nar- 
row cracks without the waist. 

The relaxation times in Eqs. ([2"5]) and ([3"0"]) do not de- 
pend on the elastic modulus of the solid matrix. However, 
we have already mentioned that for cracks with uniform 
opening, the relaxation peak (|30D can actually be formed 
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only for very highly compressible fluids (i.e., gases rather 
than liquids), for which Kf/K < a, otherwise the com- 
pressibility of rock in the crack vicinity should dominate 
and form the relaxation peak near ujk given by Eg. (|2"2"1) . 

In contrast to this, the presence of the narrow waist 
creates a strong obstacle for the fluid leaking. Therefore, 
for a given frequency, this increases the role of the liquid 
compression, whereas the pressure p msc and the deforma- 
tion of the rock in the crack vicinity is strongly reduced. 
All this results in a very strong reduction of the relax- 
ation frequency u5 r related to the compressibility of the 
fluid (compare Eq. ([50)1 and Eq. (121))) containing the addi- 
tional small factor (h/h) l/L <C 1). Thus even for "nor- 
mal" liquids like water and oil, the frequency ui r of the 
new relaxation maximum becomes comparable and even 
lower than ojk, that is, shifts from ultrasonic frequen- 
cies to the seismoacoustic range directly relevant to the 
experiment o 11 ! 13 ' 14 ^ 5 ' 16 discussed in the Introduction. 

Compared with the mean crack opening h, the local 
opening h of the waist is much more sensitive to vari- 
ations in the mean strain in the medium. Therefore it 
is not necessary to assume the existence of unrealisti- 
cally thin cracks to ensure the same sensitivity for the 
mean opening of the crack. The conditions ensuring the 
Poiseuille type of the flow are even better fulfilled in the 
vicinity of the waist, so that we can apply the consid- 
ered in the beginning of this section aguments concern- 
ing the low- and high-frequency asymptotic behavior of 
the decrement (proportional to cj and l/w, respectively). 
This means that the dissipation due to the discussed 
mechanism can be well approximated by the frequency 
dependence for a standard relaxator: 



cussed two cases of the crack without and with the waist: 



co/oj r 



1 + {Uj/Ur) 2 



(31) 



for 



where the relaxation frequency ui r is given by Eq 
a crack with a narrow waist. 

Now, using Eq. (|16|) we will determine the pref actor 9q 
in two cases: for a crack with the narrow waist and for 
a similar in size crack with uniform opening, although in 
the latter case the formation of the relaxation peak uj r 
can be possible only if Kf/K < a, that is, for gaseous 
fluids (or for unrealistically rigid solid matrix in the case 
of "normal" liquids). Nevertheless, such an expression 
will be useful for comparison. 

The velocity profile for the Poiseuille's flow in the crack 
with uniform opening h has a parabolic form 



y{h - y) 

: h 2 : 



(32) 



where v max = — ^ \ ^ \ h 2 is the maximal velocity of the 
fluid in the flow. Similarly, for the narrow waist, where 
the losses and gradients are localized, we have 



1 dP (1 



y) = -4w n 



y{h - y) 



(33) 



where v r . 



= -^l§f l^ 2 - Ec i uations OH, d32D, (BSD 

yield for dWkin/dt the following expressions in the dis- 



OWkr 

dt 



OWkr 

dt 



"3^ 



h 



(34) 



(35) 



In the estimates below we will assume that L x ~ L z ~ L. 

Let us now relate the characteristics of the flow with 
the mean strain in the crack-containing solid. When the 
mean macroscopic strain e is created in the material, the 
crack volume correspondingly changes, which creates the 
flow of the liquid inside the crack. We recall 3 ^ that the 
crack is a soft object whose relative (compared with the 
surrounding matrix) compliance £ is approximately equal 
to the crack aspect ratio, £ ~ a ~ h/L. Therefore, the 
own strain e cr of the crack (i.e., the relative variation 
in its volume) can be estimated as e cr we/aw sL/h, 
where L is the characteristic size of the crack, so that 
its volume can be estimated as V cr ~ L 2 h. Then 
for the variation in the crack volume AV cr one obtains 
AV cr « s cr V cr ~ sL 2 h ■ (L/h) = eL 3 . Under sinusoidal 
oscillatory variation of the crack volume with frequency 
uj and amplitude AV cr , we can readily relate the rate 
LoAV cr of the crack- volume variation with the maximum 
velocity of the Poiseuille's flow of the fluid inside the 
cracks without and with the waist: 



ZujeL 2 
2 h 



iujeL 2 
2~T~ 



Substituting these expressions into Eqs. (|3"4"l) and ([33)) 
for cracks without and with the waist we obtain for the 
period-averaged amounts of the dissipated energy expres- 
sions via the amplitude of the mean strain e: 



(36) 



(37) 



dWu 
dt 



dWkin 

dt 



2„2 



h 3 



-^ 2 s 2 hX 
h 3 



(38) 



(39) 



In Eqs. (|3"5|) and (|39p we took into account that for 
sinusoidal strain with amplitude e, the period- averaged 
value of its square equals e 2 /2. We also take into ac- 
count that the energy loss during one period T = 2n/to 
is Wo = (2-k / w)(dW kin / dt) aver and the density of the 
accumulated elastic energy is W e i « Ke 2 /2 (since in the 
used approximation we assume that the elastic modu- 
lus does not significantly change due to the presence of 
the cracks). Then from Eq.([T|) we obtain the following 
asymptotic low-frequency expression for the decrement 
in the case of cracks without the waists: 



9lF = 67TT7 



L 6 
h 3 K 



n cr uj, 



(40) 
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The low-frequency Eq. (|40|) is obtained neglecting the liq- 
uid compressibility and agrees with the low-frequency re- 
sults earlier obtained for cracks without waists 2 ^. 

To obtain high-frequency expressions (for ui >> u> r ) 
we take into account that because of the finite relaxation 
time, the volume of the replaced liquid is ui/ui r smaller 
than in the low-frequency limit. Thus the high-frequency 
expression for the cracks without the waists is 



Qurj 



UJ r TT ,i(fW r 

— — n cr ui — = -n cr L A ^- — . (41) 
h A K ui z 2 aK ui 



For cracks with the waists in a similar way we obtain: 



>LF 



I 



h 3 K V L 



n cr oj, 



'HF 



6nr] 



L 



I y n cr uj— 

h?K \L ui 



l / 1 



IT ^ 

2 Ucr aK u> 



(42) 



(43) 



Comparing Eqs.(|30]). and (|32"]). (g3) with the low- 
and high-frequency asymptotics of Eqs. (|3ip we obtain the 
following relaxator-like expressions for the viscous loss at 
cracks without and with the waists, respectively: 



^f 1 



uj jui r 



2 Ka 



itKjL 3 
2 Ka 



(us/uj r y 



UljuJ r 
1 + (uj/ul r Y 



(44) 



(45) 



where ui r and uj r are given by Eqs. (|29|) . 

A striking feature of Eqs. (|44|) and (|45|) is that they 
have the same maximum values determined by the char- 
acteristic size L of the whole crack. However, the relax- 
ation frequency in Eq. (|45|) is strongly shifted from ultra- 
sonic to seismoacoustic frequencies down to 10° — 10 2 Hz 
which were used in field experiment o 11 ' 15 ? 16 . 

The next point is that much smaller mean strains (such 
as the above discussed tidal strains) can already notice- 
ably affect the opening of the waist. This means that 
due to the variation in the relaxation frequency ui r for 
the crack with the waist, the position of the relaxation 
maximum can noticeably be changed, whereas the height 
of the relaxation maximum should remain yet practi- 
cally unperturbed. Taking into account the equality of 
the absolute variations Ah « Ah and the relationship 
e cr = e/a = eL/h (which we have already used above) 
we find that the relative variation in the relaxation fre- 
quency uj r caused by the variation Ae in the mean strain: 

. — i — . dui r 1 . 3/i /T .., . . 2>h 1 , . 
Aui r /uj r = As— -— Ae^(L/h) = Ae^ — . (46) 
de u) r h ha 

The parameter (3h/h)/a can be very large, for example, 
(3h/h)/a ~ 10 6 for quite realistic a ~ 10~ 4 and (h/h) ~ 
20.. 30. Thus the tidal strains with amplitude £o ~ 10~ 8 
can cause the peak-to-peak variation in the relaxation 
frequency 2Aw r /w r of several percents and, consequently, 



comparable variation 2A9/6 of the decrement should be 
observed at the wings of the relaxation curve. 

Besides, there are known indications that for flows in 
very narrow gaps (down to nanometer scale), the effec- 
tive viscosity of the liquid can significantly exceed the 
viscosity for macroscopic gap o 34 i 35 . This effect can ad- 
ditionally enhance the variations in the dissipation due 
to variations in the waist opening. Since small varia- 
tions in the average strain practically do not yet affect 
the average aspect ratio for the crack, the prefactor in 
Eq. (|4"5|) remains practically unchanged. Thus the shift 
ofJ;he relaxation maximum can cause the variation in 
AO of opposite signs depending on the position of the 
observation frequency uj relative to the frequency ui r of 
the relaxation-curve maximum similar to the situation 
for the thermoelastic loss illustrated in Fig.|3J 

Since not all cracks can have the strip- like waists, it is a 
key question how many such cracks are required to ensure 
near the characteristic frequency ui r the dissipation ~ 
10~ 2 — 10 _1 . Since the maximum of the frequency factor 
in Eq. (|4"5j) is 1 /2 we have to estimate the factor 



irK f L 3 
4 Ka 



TT K 



Tti?> (47) 



where we singled out the quantity e = L 3 n cr . The lat- 
ter is close to the effective volume of cracks 3 ^ (i.e., the 
volume of the circumscribed spheres independent of the 
cracks' aspect ratios). For the further estimates we will 
use the well known fact that the presence of cracks with 
the effective volume e results is the reduction of the elas- 
tic moduli of the material by a fraction of ae (where fac- 
tor a ~ 1 slightly differ for particular moduli 3 ^). Then 
taking for filling water Kf — 2.25 • 10 9 Pa, modulus 
K = 3.8 ■ 10 10 Pa typical of quartz, and a = 10~ 3 - 10~ 4 
we obtain 9 max « (50 — 500)?. This means that for the ef- 
fective crack density e= 10~ 2 — 10 -3 , which reduces the 
elastic moduli also by a fraction of order 10~ 2 — 10~ 3 , 
we can already obtain 8 ~ 10~ 2 — 10 _1 in the vicinity of 
the relaxation frequency uj r . On the other hand, it is well 
known that in real sandstones the modulus reduction due 
to high-compliant crack-type porosity typically is on the 
order of tens of percents 3 ^. This means that even if only 
a small portion of all cracks (a few percents or even less) 
has wavy surfaces creating the narrow waists, this portion 
can already be sufficient to explain both the background 
value of dissipation observed in the seismo- acoustic fre- 
quency range 10° — 10 3 Hz and its extremely high strain 
sensitivity indicated by the experimental dat a n i 15 ' 16 . 

Therefore, the considered modified mechanism of the 
squirt-type dissipation in cracks suggests a plausible 
alternative explanation to the experimentally observed 
rather high dissipation in the seismoacoustic frequency 
range. In terms of paper 3 -, our mechanism is of purely 
"microscopic" type and does not require the presence 
of larger-scale (" mesoscale" ) heterogeneities assumed in 
the "double porosity" and "patchy saturation" models 
to shift the frequencies of viscous relaxation towards the 
seismoacoustic frequency range. Certainly such mecha- 
nisms can operate simultaneously. However, among those 
possibilities, only the above considered modified squirt 
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FIG. 5. Schematically shown relative positions and heights 
of the viscous relaxation peaks at LO r and ujk for the crack as 
a whole (curves 1 and 1 ') and the low-frequency peak for a 
crack of the same size L having an inner narrow waist with 
the local opening h and length I (curve 2). The examples 
correspond to (h/hf(T/L) = 10 4 and a = 10 3 . 



strong peaks in the seismo-acoustic frequency range, so 
that even for small density of such cracks, their contribu- 
tion can easily account for the typically observed levels 
of the dissipation in this frequency range. 

Concerning the question of averaging over the distribu- 
tion of real cracks over their parameters, we can put for- 
ward very similar arguments as for Eq. in the above 
considered case of thermoelastic loss. Namely, it is rea- 
sonable to assume that the distributions over the sizes L 
and h of the crack as a whole and the distribution over the 
local parameters h and I of the waist (which determine 
the characteristic relaxation frequency ui r ) are essentially 
independent and can be factorized. Therefore, the aver- 
aging over the characteristic size L of the crack gives 
only a numerical factor like in integral Eq. (jlip , whereas 
the averaging over the relaxation frequencies (of, equiv- 
alently, the relaxation times) can also be performed in- 
dependently. Such averaging should not radically change 
the conclusions obtained for the simplest case of cracks 
with identical parameters of the inner contacts or waists. 



V. CONCLUSION 



mechanism is evidently able to ensure sufficiently high 
strain sensitivity to explain the observation o 1 1 ' 1 5 ' 16 . 

To better understand the ratio of contributions of the 
conventional and the proposed modified squirt mech- 
anisms, let us compare the corresponding decrements 
given by Eqs. (|4"4"|) and (|4"5)) in the vicinity of the char- 
acteristic relaxation frequency u) ~ ui r of the cracks with 
the narrow waists. Since uj r <C ui r for this comparison we 
take the low-frequency asymptotic form of Eq. (I44|) , such 
that in the vicinity of ui r we obtain 

0( w ~5 r )«2^^0 ma:c , (48) 

where the difference in the relaxation frequencies U) /ui r — 
(h/h) 3 {T/L) »1 as follows from Eqs. (J23J and (J3DJ). The 
factor (h/h) 3 (l/L) which can easily reach 10 2 — 10 4 , so 
that quite a small portion of cracks with narrow waists 
can ensure in the low-frequency range a contribution 
comparable with or even strongly exceeding the con- 
tribution of the majority of other cracks of the same 
size, but without the waists. For example, a portion 
n cr /n cr ~ 1CP 2 — 10~ 4 of cracks with narrow waists is 
able to ensure extremely high sensitivity of the overall 
decrement to weak variations in the mean strain. 

Figure [5] illustrates the relative positions and heights 
of the relaxation peaks corresponding to the relaxation 
frequencies u) r (Eq. (p^ ) for a crack with a narrow 
waist, and the peaks at frequencies u>k (Eq. (|22|l) and u> r 
(Eq. for cracks with the same aspect ratio a and size 
L, but without waists. We emphasize that for lok < w r , 
the relaxation peak at uj r actually does not exist and 
is shown by the dashed line. However, it is shown for 
convenience of comparison with the low-frequency peak 
u5 r , for which the height is the same as for the would-be 
peak at u> r . The figure demonstrates that the viscous 
relaxation in cracks with the waists can form extremely 



The performed analysis of the role of elongated inner 
contacts and waists in cracks significantly changes the 
conclusions based on conventionally discussed models of 
thermoelastic dissipation at cracks (lik o 23 ' 24 ) and viscous 
squirt loss (like work a 28 ' 29 ' 31 ' 33 ' 38 ). Thus a single larger 
crack with a strip-like contact can ensure the same ther- 
moelastic dissipation as 10 5 — 10 6 small cracks of the size 
equal to the contact width. 

For fluid saturated cracks, a rather intense maximum 
formed by a crack with a narrow waist can ensure the 
same dissipation in the seismo-acoustic frequency range 
of 10 2 — 10 3 Hz as a similar crack without the waist would 
produce in the ultrasonic range according to conventional 
squirt-dissipation models (see Fig. [5]). This can signifi- 
cantly affect some conclusions^ 8 on insignificant role of 
local viscous loss at cracks in the seismo-acoustic range. 

Probably the most striking feature of the considered 
modified dissipation mechanisms is their giant strain sen- 
sitivity. In this context, it should be clearly understood 
that each group of cracks with the wavy asperities can ex- 
hibit this giant strain sensitivity only in a rather narrow 
strain range. When the waist becomes either completely 
closed or widely open, the loss at such cracks does not 
much differ from that at cracks without the asperities. 
Nevertheless, since the parameters of real cracks should 
have rather wide distribution, for a current mean strain, 
another portion of cracks with such narrow waists or con- 
tacts can be "activated". This resolves the problem 14 
how apparently very soft defects exhibiting giant strain 
sensitivity can exist under very different mean pressures 
(in a wide range of depths in field conditions). Their gi- 
ant effective compliance should be understood as differ- 
ential. It cannot be directly extrapolated for significantly 
higher strains 10~ 6 — 10~ 4 often used in laboratory exper- 
iments or typical of tectonic strains, for which the level 
of the dissipation does not change many times, although 
for 10 -8 the variations may reach several percents. 
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